cd 



Peak positions and shapes in neutron pair correlation functions 
from powders of highly anisotropic crystals 



O 
O 

o 

(N 

j^' Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

C/3 



D. A. Dimitrov, H. Roder, and A. R. Bishop 

laboratory, Los Alan 
(February 1, 2008) 



O 

&> ; Abstract 



The effect of the powder average on the peak shapes and positions in neu- 
tron pair distribution functions of polycrystalline materials is examined. It is 
shown that for highly anisotropic crystals, the powder average leads to shifts 
O ' in peak positions and to non-Gaussian peak shapes. The peak shifts can be 

as large as several percent of the lattice spacing. 
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I. INTRODUCTION 

The purpose of this paper is to demonstrate that the peak shapes in the powder averaged 
static pair distribution functions (PDF) of highly anisotropic harmonic crystals are generally 
non-Gaussian and the peak positions (maxima) are at different distances than the average 
atomic pair distances. Our approach is to consider the effect of the powder average on 
specific peaks of the theoretical PDF. We discuss the relevance of these results to what is 
actually measured in neutron diffraction experiments. 

The lattice dynamics properties of highly anisotropic crystals, in particular for crystalline 
materials of layered and chain structure, are known to lead to significant deviations from 
the Debye approximation for thermodynamic quantities rendering it inapplicable in these 
cases. LifshitzEl showed that there are different low temperature regimes in which the specific 
heat does not behave as T 3 but rather as T a with the exponent a having different values 
depending on the different regimes for the phonon dispersion relations. The behavior of 
the thermal expansion tensor is similar. The anisotropy also leads to very different in 
magnitude Debye- Waller factors^ and to violation of the Cauchy relations for the elastic 
constants. The latter is due to the inability to describe the inter-atomic forces with central 
pair potentials. A number of different models have been developed for the inter-atomic 
potentials to deal with the anisotropy, depending on the physical properties of the crystal 
concerned. For example, generalized force constants models are considered for some simple 
metals, bond bending forces are introduced for the strongly anisotropic covalent crystals, 
anisotropic ionic polarizabilities for ionic crystals, etc. It is of interest to determine the 
effects of highly anisotropic harmonic lattice vibrations together with the powder average 
on the static pair distribution function (PDF) and the implications for the corresponding 
measurements in neutron powder diffraction experiments. 

The space-time pair correlation function G (r, t) introduced by Van HoveS in the theory 
of neutron scattering has been primarily studied, both experimentally and theoretically, 
for liquids and amorphous materials.^ The correlation function in liquids is often assumed 



to be Gaussian provided the system is considered as isotropic.™! The interpretation of the 
experimental data is quite simplified by assuming Gaussian peaks, and indeed much has been 
learnt about the structure and the interactions in such systems from the pair correlation 
function. 

The determination of the static PDF from experimental data on powder samples is ob- 
tained by first extracting the static structure factor S(q), where q = |q|, from the measured 
(effective) coherent cross section ( ^'' ) ff and then Fourier transforming it. S(q) has struc- 
ture for large momentum transfers hq in crystals due to their long-range periodic structure. 
For an accurate PDF to be obtained from the Fourier transform of S(q) in such materials, 
the structure factor should be known to large enough values of q. Only recently the de- 
velopment of high-intensity neutron sources has allowed the measurement of effective cross 
sections up to q ~ 50 A" 1 or larger and thus the ability to obtain a PDF for polycrystalline 
materials from neutron diffraction experiments.Q'B 

The extraction of the structure factor S(q) from the neutron scattering data is based 
on several assumptions. The first Born approximation is assumed to be applicable for the 
interpretation of the differential cross section. This is the case, at least in the static ap- 
proximation when the energy transfer of the scattered particle or photon, in the case of 
x-rays, is negligible compared to its own energyB The static approximation is applicable 
to the scattering of x-rays but cannot be applied directly to the scattering of neutrons be- 
cause of their specific energy and wavelength scales (note that the static structure factor 
S (q) which is obtained by integrating out the u dependence of S (q, u>) is different from 
the structure factor obtained in the static approximation!). Corrections to the static ap- 
proximation have been developed by PlaczeklI9 and by WickO which are applicable for high 
incident neutron energies. Yarnell et alx3 have shown that the application of the Placzek 
correction to neutron detectors with different efficiency can lead to determination of the 
structure factor S(q) with an accuracy of ~ 0.01 for liquid Ar. However, it has been pointed 
out by Ascarelli and Cagliatil^l that there are cases when the Placzek correction could lead 
to large errors, particularly for values of q for which the curvature of S(q) is high. Generally, 



this requires a careful study of the Placzek correction before its application is adopted. We 
shall assume that the Placzek correction is applicable in order to determine the S(q) from 
the differential scattering cross section, as is usually done when analyzing the experimental 
data&3 for liquids as well as for crystals.^ We also assume that all the other corrections and 
normalizations applied to the raw data, such as multiple scattering, incoherent scattering 
calibration, absorption, polarization, etc., and measurements to large momentum transfers 
are made so that the Fourier transformation termination errors are very small. A careful 
experimental determination of the PDF which satisfies these assumptions will allow for a 
meaningful comparison with a theoretical calculation.^ 

In Sec. U, we introduce the static pair distribution function for harmonic crystals with 
arbitrary structure and its relation to the experimental quantities measured. The powder 
average for specific cases is considered in Sec. |I fl| which contains the main results of this 
paper. The implications of these results together with their relevance to x-ray experiments 
are discussed in Sec. [TV]. 

II. EXPERIMENTAL AND THEORETICAL BACKGROUND 

The coherent differential cross-section for scattering of thermal neutrons in the first Born 
approximation is given byO 
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where a co h = ^iib (b is the compositionally averaged scattering length!^), k and k' = k — q 
are the initial and final wave vectors of the scattered neutron, r^ (t)'s are the operators for 
the position vectors of the iV particles in the target system in the Heisenberg representation, 
S (q, uS) is the dynamical structure factor, and (...) denotes thermal average. We have here 



assumed that the scattering system consist of only one type of atoms. In the general case, 
the different types of atoms will have different scattering lengths and this should modify 
Eqs. (|IJ) and (|2|) accordingly. This case will be considered explicitly for the static PDF in 
crystals. 

The Van Hoveta pair distribution function is related to S (q, to) via a space-time Fourier 
transformation 

G(r,t) = -±xjs (q, w) e-^^dqdt (3) 

This function is complex since in the general case the position vector operators for the 
particles of the scattering system do not commute at different times. At equal times, t = 0, 
these operators commute and G (r, 0) is real. 

Provided S"(q, uj) can be extracted from the measured differential cross-section, Eq.([I]), 
then Eq. (|IJ) can be used to obtain G (r, t). The pair distribution function has a simple 
physical meaning and is easier to understand than S (q, uj). It is not currently feasible 
to determine experimentally the dynamical structure factor in such a domain of the four 
dimensional (q, uj), space for the Fourier transformation in Eq. (^) to be applied. However, 
in some cases the uj integration can be effectively done by the neutron detectors. This allows 
for the determination of the static structure factor 
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S (q, uj) duj = G (r, 0) e iq - r dr, (4) 
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provided q is kept constant in the u integration. This is the case, at least, in the static 
approximation. The neutron detectors carry out this integration at constant scattering 
angle 9 and the effective coherent cross-section per unit solid angle is 

\dn) coh te 2rfi./_ 0o ^ >k y 

where hw max is the energy of the incoming neutron, e(k') is the detector's energy dependent 
efficiency, and q is a function of the scattering angle 9 via q 2 = k 2 + k' — 2kk' cos(2#). The 
problem is to obtain 5(q), as defined by Eq. (|J), from the measured effective cross section 



(m) h w ^ en the static approximation is not valid. The usual way to tackle this problem 
is to apply the Placzek correction as mentioned in Sec. [|, particularly when the incident 
neutron energy is high compared to the energy transfers in the scattering processes. 

For powder samples, the correlation functions do not depend on the angular coordinates 
of q and r since they are effectively averaged over these variables. Then the relations between 
S(q) and 

p( r ) = G(r, 0) - 8(r) (6) 
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where po = N/V, V is the volume of the system, and the forward scattering has been 

subtracted from S(q), which is the conventional practice since it is negligible for scattering 

angles larger than ~ 10~ 2 seconds of arc.t3 The experimental PDF, p exp (r), is obtained from 

Eq. (|S|) with the upper limit of the integration q = q m ax determined from the maximum 

value of the scattering vector accessible in the measurement. If the assumptions made in 

Sec. | regarding p exp (r) are valid, we can compare this with a theoretically derived p(r) 

which we consider next for harmonic crystals. 

The static pair distribution function, p(r) , in a form suitable for neutron diffraction from 

a crystal can be defined by 
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where Nf, is the number of basis atoms in the unit cell, R(/) is the position vector of the 
origin of the Ith unit cell, d, d' are basis atom indices, b(ld) is the scattering length of the 
atom with equilibrium position at H(ld) = R(Z) + R(c/). H(d) is the position vector of 
the dth atom in the basis relative to the origin of the cell. The r(0d; I'd') = r(0d) — r(l'd') 



is the difference between the instantaneous positions of the two atoms. The instantaneous 
position of an atom is r(ld) = R(W) + u(ld), where u(ld) is the deviation vector from the 
equilibrium position (the time arguments of the two u vectors in Eq.(^) are not explicitly 
written since they are at equal times). The sum excludes the terms with /' = 0, and d = d! . 
For harmonic crystals 
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where R(0d; I'd') = R(0d) - R(Z'd'), and the tensor T(0d; I'd') is 

T af3 (0d; I'd') = /(u (Od) - u (l'd')) a (u (Od) - u (Z'd'))^} • (H) 

The result for G (r, t), t =£ 0, in such crystals is given by Van Hove. til The thermodynamic 
averages of the type (u Q (Id) up (I'd')) are easily evaluated once the phonon frequencies uj s (k) 
and polarization vectors cr s d (k), where s is the branch index, are known for reciprocal space 
k-vectors in the first Brillouin zone 
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Here n s (k) = (exp (/3huj s (k)) — 1)~ , N k is the number of unit cells in a crystal with periodic 
boundary conditions (or k points in the first Brillouin zone), and Ma is the mass of the 
dth type of atom. Note that (u a (Id) Up (Id)) is related to the Debye-Wallerc3 factor via 
2H / rf(q) = q • (u (Id) u (Id)) ■ q, where q is a reciprocal space vector. 

As expected, the static PDF for harmonic crystals as a function of the vector argument 
r, Eq. (|I0"D, consists of a sum of Gaussian distributions for each pair of atoms, apart from 
the scattering length factors. Each distribution is centered at the pair's average distance 
R(0d; I'd') and its second moments with respect to R(0d; I'd') are given by the elements of 
the tensor T(0d; I'd'). These moments determine the width of the peaks in the PDF along 



any direction in r-space. It is clear from Eq. (|12|) that the moments depend on R (/) — R (/') 
and the type of atoms in the pair. Since (u a (Id) Ug (I'd')) — * when |R (/) — R (l')\ —>■ oo, 
the moments of the distribution will approach (u Q (Od) Up (Od)) + (u a (Od') up (Od')) and will 
become distance independent in this limit. 

III. POWDER AVERAGE 

When the scattering is from a crystalline powder sample, the structure factor depends 
only on the magnitude of the scattering vector and is averaged over its angular coordinates. 
Then the pair distribution function defined by Eq. (|j), or by the inverse Fourier transform 
of Eq. (f|) for the static PDF, will depend only on r and not on the direction of r. In the 
latter case, this is equivalent to an angular average in real space 

P(r) = f f^pW, (13) 

JSn r 4vr 

where r = |r| and Sci r is the surface of a sphere with a center at the origin of the coordinate 
system and a radius r. The static, angularly averaged, S(q) and p(r) are related by Eqs. 
(0) and (|). 

Below, we will apply Eq. ( |13"D to the different peaks in p(r), as given for crystals by 
Eq. (pi]), in order to examine the effects of the powder average on the static PDF. Note 
that if p(r) is to be obtained from S(q) via Eq. ([§]), the integral over q must be carried out 
to infinity rather than to a maximum value q ma x determined by the scattering instrument. 
In what follows, it is undedstood that the assumptions made in Sec. | hold such that the 
theoretical PDF obtained from Eqs. ( |10D and ([TB]) can be compared with the experimentally 
determined. 

We consider first the case of an isotropic tensor T(0d; I'd') because it is an important 
limiting case against which we can compare the case of a general form of T. 



A. Gaussian Behavior of rp(r) Peaks 

The simplest case for which the powder average integral, Eq. flTS|) , can be taken analyt- 
ically is for an atomic pair with an isotropic T(0d; I'd') = t 2 X tensor (note that generally 
in taking the integral in Eq. (|l3|) for each pair, the r-space coordinate system can always 
be rotated such that T(0d;l'd') = diag{tj, t 2 ,, t|}), where X is the identity tensor. The 
contribution P(od-,i'd')( r ) of such a pair to the total PDF is 

b(0d)b(l'd') 
P( ° w(r) = lHnrR^d ] l'd') {PNWm ' d ' ) ' t2) (r) 

-PN(-R(Od;l'd'),fi)(r)), (14) 

where PN(R,t 2 ) ( r ) — e~^ r ~ R ' ' 2t /V2nt 2 is the density of the normal distribution with a mean 
equal to R and standard deviation t. The contribution of PN(-R(od;i'd'),t 2 ){ r ) to the PDF can 
often be neglected when t 2 is large enough, and for r of the order of the pair distances in 
a crystalline solid which are usually equal to several Angstroms or larger. Then the peak 
in rp(r) due to such a pair of atoms will be approximately a Gaussian (centered at the 
equilibrium distance R(0d; I'd') between the atoms and with a standard deviation equal to 
t) divided by R(0d; I'd') (apart from the constant scattering length factor b(0d)b(l 'd') / 'Anb 2 ) . 
Therefore, if for all pairs T (Od; I'd') = t 2 (Od; I'd') I were isotropic (which is not the case for 
crystalline solids) the powder averaged PDF would be 

/ s J_ v^ , b{Qd)b{l' d') PN{R{od;Vd'),t 2 {od-Vd')){r) 
n(i>),d,d> v ' ' 

Strictly, 4irrp(r) is not a sum of Gaussians (again neglecting the scattering length factors) 

since ea.ch.p N ( R ( d t i'd'),t 2 {od;i'd'))( r ) is also scaled by its mean value R(0d; I'd') which is different 
for the pairs nonequivalent by symmetry. The form for an isotropic PDF peak given by 
Eq.(|T4"D can be used to approximate peaks in diagonally cubic crystalsE-3 when 

\(u a (Od)u (l'd')) + (u a (l'd')u p (0d))\ < u 2 {0d) +u 2 {l'd'), 

for a, (5 — x, y, z, since in such crystals (u(ld)u(ld)) = u 2 (ld)T. This may be applicable 
to pairs of atoms separated by a large distance relative to the nearest neighbor distance 

9 



such that the correlation terms become much smaller than the diagonal elements of the 
self correlation ones. This is a special case of the more general statement that the long- 
range part of the PDF can be obtained to a very good approximation only from the Fourier 
transformation of the elastic part of S(q). 

B. Non-Gaussian Behavior of rp(r) Peaks 

We consider next a case which allows us to study the effect of the powder average for 
a strongly anisotropic tensor T and which can still be handled analytically. The powder 
average contribution to the PDF due to a pair of atoms with T(0<i; I'd') = diag < t±, t\, t 2 >, 
t\ < t 2 ± , and R(0d; I'd') = (0, 0, R) is 



PW#) (0 = 4^: {PN(R,tl) ( r ) G ^ ~ hR ) + P N(~R,tl) ( r ) G ( flr + hR ) ) > ( 16 ) 

where 



G(x) = Jirt 2 /2t 2 ± (t\ - tfj exp (x 2 ) erf (x) 



a = \j it]_ — t\j /2£j_£m, b = \ t 2 L /2t 2 (t 2 ± — t 2 ,), erf(x) is the error function, and we have 
omitted the scattering lengths scaling factor b(0d)b(l'd')/b 2 . It is straightforward to check 
that in the limit t 2 —>■ t\ from below we recover the result for the isotropic case given 
by Eq. (|14]) . These specific forms of T and R arise, for example, for specific 0\-B or 0\-0\ 
nearest neighbor pairs of atoms in ABO3 cubic perovskites. We use here Cowley's notation^ 
for this structure. The self correlation term (u(B)u(B)) for the B atoms is isotropic and 
the anisotropy in the T is due to the oxygen self correlation term. The latter is proportional 
to its Debye- Waller factor and has the same form as T. The B-0 correlation terms also 
contribute to increasing the anisotropy of T. For the 0\-0\ pair of atoms, the T tensor is 
even more anisotropic because of the contributions of both atom's anisotropic (u(Oi)u(Oi)) 
terms. Physically, this is associated with the local environment of the 0\ atom. It has four 
A nearest neighbors (nn) in the xy plane and two nn B atoms along the z axis. The oxygen 

10 



polarizability is also very anisotropic and it depends strongly on the oxygen's environment xB 
To study the peak positions and shapes as functions of t]_ and ti we have chosen R = 1.956 
and R = 3.912 A(the lattice constant in cubic ABO3 perovskites is often close to 3.9 A) which 
corresponds to the first nearest neighbor 0\-B and 0\-0\ bond distances. The t\, t?, are 
varied in intervals such that 5 < t^ 2 < 99 A~ 2 and 100 < £7 2 < 250 A~ 2 . Lattice dynamics 
calculations^ on Lao.7Sro.3Mn03 using the cubic perovskite structure and the nonlinear shell 
model of Migoni et alx3 show that these intervals for t 2 ± , t 1 could be physically relevant. The 
peaks of the normalized (to unity) distributions p(r) = Const x rpn{r) for two ( R, £]_ , £JT 
sets of parameters, (1.956, 21.71, 133.99) and (3.912, 6.23, 248.11), are shown in Fig. |. 
These values of of (R, t^ 2 , t~^ 2 ) arise for a specific set of the shell model parameters chosen 
from a Monte Carlo sample sequence in a domain of their phase space. The sequence was 
generated in a Reverse Monte Carlo estimationll3 of the shell model parameters from PDF 
data.§ 
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FIGURES 
FIG. 1. Non Gaussian peaks in the normalized distribution of rpn(r) vs r/R for two values 

of ( R, £7 , t~ ) • The magnitures of R are given in the inset of the figure. Each peak position is 



shifted at a value greater than the corresponding average distance R. 
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We calculated these peaks using a numerical implementation of Eq. ( |16D with a step 
in r equal to 10~ 4 A. The maximum of each peak (its position) is shifted to a distance 
larger than the corresponding average R. This shift is equal to 0.0186 A for the peak with 
R = 1.956 A and to 0.0321 A for the one with R = 3.912 A. These are fairly large peak 
shifts and, importantly, of the same magnitude as observed peak shifts in specific materials 
interpreted as statitic displacements from the periodic lattice structure. The peak shapes 
show deviation from the Gaussian, as the asymmetry in the shape with respect to the peak 
position is clear from the figure (a quantitative measure for this will be given below). The 
p(r) decreases to zero slower for r > R ma x, where R ma x is the distance of the peak position, 
than for r < R ma x- 

It is of interest to note a recent studycll on the PDF of the superconducting La 2 _ x Sr a; Cu04 
in which a change of the Cu-0 bond length equal to 0.024 A was reported. However, this 
bond change has been attributed to a different mechanism. This is a layered material and 
it will clearly be of importance to be able to estimate to what extent this shift may actually 
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be due to the effect of the powder average. 

We proceed to formally study the peak shift relative to the average R and the deviation 
of the peak shape from the Gaussian for the domain of (£J_ 2 , £7 2 ) given above and for 
R = 1.956 A. The peak shift 
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R 



x 100, 



(17) 



change relative to the average distance between the atoms in the pair is plotted as a function 
of t~ 2 and t|T 2 in Fig ||. 

FIG. 2. Contour plot of the peak shift magnitude, Eq. (|17|) as a function of the £7 and t~7 



parameters. 
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The figure shows that increasing the anisotropy of the eigenvalues of the T tensor could 
lead to peak shifts as large as 3% of R. Moreover, the contour plot gives that for shifts larger 
than approximately 0.5%, the value of the £7_ 2 remains almost unchanged on each contour 
level while increasing t^ 2 . Note that peak shifts of the order of 0.5% (~ 0.01 A for R = 1.97 
A) or larger should not be difficult to detect experimentally. In the limit t~ 2 — > t^ 2 from 
below, the peak shift tends to zero as it should in the isotropic limit and the peak position 
is at the average distance R. 

The deviation of a peak in rp{r) from a Gaussian form can be studied^ by calculating 
the moments of the normalized (to unity) distribution given by this peak. If the normalized 
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distribution were a Gaussian p7v(r,o- 2 ), the moments obey 



(r — r) = 0, 



(r-r) 2n = (2n-l)\\((r-r) 2 \ 



for n — 1, 2, . . ..When the normalized distribution deviates from a Gaussian, the quantity 

(18) 



o. 



(r — r) 



-\2n 



(2n-l)H (r-r) 



n 1, ^ — A >Jj • • • j 



will deviate from zero. The contour plots of the dependence of a n (n = 2, 3, 4) on ( £,, 2 , t ± 2 j 
for the same intervals as in Fig. ^ and for the same peak is given in Fig. ^. 

FIG. 3. Contour plots of the dependence of a n , (n = 2, 3, 4) (see Eq.(^q)) on f tjT , t~]_ ) for 
the same peak as in Fig. g. 
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The behavior of a n (n = 2, 3, 4) as a function of ( tj 2 , t^ 2 ) is very similar to the behavior 
of the peak shift shown in Fig. 0, as expected. When T is close to isotropic the peak shift 
is less than 0.1%, Figs. and |3|, which is < 0.002 A for this peak, and the peak shape is 
close to a Gaussian (a n < 10~ 5 , n — 2, 3, 4). The values of «2, «3, and a^ are all positive 
for these intervals of t^ 2 and £J_ 2 indicating that the normalized distribution of rp{r) for this 
peak tends to zero slower than a Gaussian with increasing r. 

The magnitude of the peak shift and the behavior of the a n quantities plotted in Fig. 
H] and Fig. [| depend on the anisotropy of the tensor T, e.g. determined by the ratio of 
its largest to smallest eigenvalues, and on the average pair distance R. Fig. || shows that 
for t~ 2 < 40 A -2 , the peak position shift is approximately constant for a constant £7_ 2 and 
varying t^ 2 . Thus, the peak shift of 0.82 % for the peak with R = 1.956 A plotted in 
Fig. [I] and a ratio of t^ 2 /t^ 2 fh 6.17 could be obtained for the same R and t|T 2 /tJ 2 « 4. 
Such ratios, and even much larger, have been experimentally observed at least for the self- 
correlation tensors (uu), as can be seen from published Debye- Waller factors in highly 
anisotropic materials. For example, ( u 2 ) / (u\) ~ 3 has been found for layered graphite 
crystalsE3 and (u 2 ) / (u 2 ) ~ 19 for the moniclinic 7M structure! 2 ^ in Ni62.5Al37.5- 



The peak shift for a peak given by Eq. fllq ) decreases when increasing R. This 
shift was 0.0186 A for the peak due to the 0\-B nn pair of atoms, (R, £7_ 2 , £7 2 ) = 
(1.956, 21.71, 133.99), and it reaches 0.001 A fori? ^37 A and the same values of (ij_ 2 , t7 2 
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While the eigenvalues of the T tensor will also change with R due to the pair displace- 
ment correlation contributions to T, these contributions will decrease with increasing R and 
for large distances, the eigenvalues of T will be approximately given by the displacement- 
displacement self correlations. Then for very anisotropic Debye- Waller factors T will also be 
very anisotropic and peak shifts will result for the powder average. More importantly, these 
shifts will be different for the atomic pairs at different distances with shifts as large as e.g. 
0.02 A for the nn bonds and decreasing slowly to 0.001 A for R rs 37 A. Moreover, these 
shifts are confined to the pairs in which at least one of the atoms is O, in our perovskite 
example, while peaks due to A-A or B-B pairs could be positioned very close to their respec- 
tive average i?'s at all distances, because for these atoms the Debye- Waller factors are very 
close to isotropic. Thus, we may have a crystal with several different types of atoms in the 
unit cell, very well defined periodic structure, and nevertheless some of the peak positions 
in its powder averaged PDF could deviate one or more percent from the average atomic pair 
distances while other peaks could be practically positioned, within the resolution in r space, 
at the average i?'s. 

IV. SUMMARY 

We have considered the effect of the powder average on the peak shapes and positions 
of the pair distribution function as measured from neutron diffraction, assuming that the 
harmonic approximation is a valid description for the phonons in the system. When the pair 
correlation displacement tensor T(ld; I'd') for a given pair of atoms is isotropic the shape 
of the peak in rp(r) (and not in p(r)) corresponding to this pair is a Gaussian scaled by 
the average distance R(ld; I'd') between the atoms and a factor which depends on the scat- 
tering lengths. However, for polycrystalline samples, the pair correlation tensor T(ld; I'd') 
is generally not isotropic. Then the powder average leads to deviations from the Gaussian 
shape and to peak shifts depending on the "degree of anisotropy" of T and the magnitude 
of the average distance R. The anisotropy of T depends on the relative magnitude of its 



16 



eigenvalues. For the specific T and R of a pair of atoms, the magnitude of the peak shift can 
be as large as several percent of the of the magnitude of the average pair distance R. This 
could be of the order of 0.02 A or larger. Such a change in a peak position should be possible 
to detect experimentally given the currently available, high resolution neutron diffraction 
measurements using high intensity neutron pulse sources. Analysis of the moments of the 
normalized distribution of a peak in rp(r) on the anisotropy of T shows that the normalized 
distribution tends to zero slower than a Gaussian for increasing r. It is also possible to have 
a highly anisotropic material for which some of the peaks in its powder averaged PDF are 
displaced markedly from the average positions while other peaks are at their atomic pairs 
average distances. In this case, the Rietveld analysis of the Bragg peaks may show that the 
crystal is periodic. 

Generally, for any experiment which allows to obtain the powder averaged p(r) accurately 
enough for (perfect) harmonic crystals in order to make a meaningful comparison with the 
theoretical PDF as obtained from Eqs. ( |T0| ) and (|T3|), the peak shifts and shapes reported 



here should be relevant effects. We considered in this paper the case of neutron powder 
diffraction experiments which we expect to be an appropriate probe to detect the peak shape 
and position effects we predict based on the theoretical results. Other experiments which 
allow to probe the PDF of a poly crystalline material are x-ray scattering and EXAFSHj 
Similarly to the neutron experiments, a structure factor can be deduced from x-ray scattering 
data. In principle, from its inversion (Eq. (||)), a PDF can be obtained. This is, however, 
more complicated for x-rays than for neutrons since the scattering lengths for x-rays (the 
atomic form factors) are q-dependent rather than just numbers in the case of neutrons. 
Powder averaged PDF are often used in the interpretation of EXAFS spectra but Gaussian 
peaks are assumed at the average distances. 

Clearly, experimental confirmation of our result will be very important for the proper 
analysis and understanding of PDF data in very anisotropic materials. 

Work at Los Alamos is supported by the Department of Energy under contract W-7405- 
ENG-36. 
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